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We study the thermodynamical properties of a mesoscopic 
Fermi gas in view of recent possibilities to trap ultracold 
atoms in a harmonic potential. We focus on the effects of 
shell closure for finite small atom numbers. The dependence 
of the chemical potential, the specific heat and the density 
distribution on particle number and temperature is obtained. 
Isotropic and anisotropic traps are compared. Possibilities of 
experimental observations are discussed. 



investigate the case of an isotropic harmonic trap. We 
analyze the influence of the shell structure on the chemi- 
cal potential and the specific heat as a function of number 
and temperature. For small particle number, density dis- 
tributions deviating from the Thomas-Fermi distribution 
are obtained. In Sec. Ill the anisotropic trap is consid- 
ered with respect to chemical potential, specific heat and 
density. The deviations from the isotropic trap are dis- 
cussed. In Sec. IV the results are summarized in view of 
experimental realizations, and the validity of the assump- 
tions underlying our theoretical approach is discussed. 
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I. INTRODUCTION 

The recent realizations of Bose-Einstein condensation 
in dilute atomic vapors have not only stimulated 

many investigations on Bose atoms but also studies of 
degenerate Fermi gases. As opposed to charged Fermi 
gases, i.e. nucleons or electrons in solids, the effects of 
Fermi statistics in neutral atomic gases occur at much 
lower temperatures (typically below 10~ 7 K) and at den- 
sities which allow a treatment as dilute quantum gas be- 
cause of the weak interatomic interactions. Due to the 
selection rules for collisions, spin-polarized fermionic al- 
kali atoms, like 6 Li or 40 K in magnetic traps, remain 
metastable in the regime of quantum statistical degener- 
acy. Due to Fermi statistics the lowest scattering channel 
(s-wave scattering) is closed for atoms in identical mag- 
netic sublevels. Therefore an ultracold spin-polarized 
Fermi gas will be less influenced by interactions than the 
ultracold Bose gases Q-^] . Also a BCS transition as stud- 
ied in j4|J^] is naturally excluded here. The only remain- 
ing interaction is the magnetic dipole-dipole interaction 
between the atoms. An estimate of its contribution to 
the mean- field § yields (V DD ) = {ha d /M)( r - 3 } where 
M is the atomic mass and ad < 10~ 10 m. This effect is 
neglected here. 

The purpose of the present paper is rather to consider 
the stationary features of an ideal Fermi gas in isotropic 
or anisotropic harmonic traps. Since our results are based 
on the numerical calculation of the state sum without 
further approximation, they are complementary to the 
recent paper of Butts and Rokhsar |Q, where a contin- 
uous spectrum and Thomas-Fermi approximation were 
used. That treatment becomes exact in the limit of large 
particle numbers. In contrast, we here focus on the ef- 
fects of small particle numbers, where the shell structure 
still affects the behavior of the many-particle system. 

The outline of the paper is as follows. In Sec. II we 



II. PROPERTIES OF THE FERMI GAS IN AN 
ISOTROPIC HARMONIC TRAP 

A. Degeneracies and the Fermi edge 

We first study the isotropic harmonic trap with poten- 
tial energy 



Mu 2 , , , 2 \ 

V = — - (x 2 +y 2 + z 2 ) 



(2.1) 



and frequency u>, because of its distinct features com- 
pared with the anisotropic trap studied below. In this 
case the degeneracy of states with equal energy 



is given by 



lu = -{v+l){u + 2), 



(2.2) 



i.e. equal to the number of simple partitions of v as a 
sum of three integers v = v x -\-v v -\-v z . Since <?„ gives the 
degeneracy of a shell of energy E v , one finds the total 
number S a of quantum states with energy smaller than 
E a as the sum over the shells < v < a, 



y^g^ = s a . 



u=0 



The sums 



S a 



6 



(a+l)(a + 2)(a + 3) 



(2.3) 



(2.4) 



define the sequence {S a } — {1, 4, fO, 20, 35, 56, ...} and so 
forth. Note that each oscillator state is assumed to be 
filled with a single Fermion since only one spin orientation 
is confined by the magnetic trap. 
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For simplicity, our calculations are done using the 
grand canonical ensemble B. The thermal occupation 
of a state with energy E v at a temperature k B T = 1//3 
is given by Fermi-Dirac statistics as 



1 



z 1 exp(PhLuiy) + 1 ' 
where the fugacity z is determined from the condition 



(2.5) 
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g v n v = N. 



(2.6) 



The definition of the fugacity z = exp/3(tt — (3/2)hui) 
absorbs the zero-point energy. N is the total number of 
particles in the trap. For a given particle number one 
can determine the Fermi energy 



E f 



(2.7) 



where vp is the shell up to which the trap levels are 
filled with particles at temperature T = 0. In this limit 
(P — > oo) the Fermi-Dirac distribution approaches a step 
function. Note that in the case of a mesoscopic ensemble 
the zero-temperature equation 
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v=0 



g v 0{v F -v)=N 



(2.8) 



does not have a solution for each N, but only for the 
discrete set of total particle numbers N £ {S a }. For 
N £ {S a } we may still define vp = \xp~\ as the small- 
est i nteg er equal or greater than the exact solution of 
Eq. (ET 



xp = A 



1 

3A 



where A is given by 

A = (3N + v/97V 2 - 1/27 



1/3 



(2.9) 



(2.10) 



The expression Ep — hu>{\xp ] +3/2) has to be compared 
to the Thomas-Fermi approximation for the Fermi edge, 



huj(6N) 1/3 . 



(2.11) 



B. Calculation of the chemical potential 

We now turn to the determination of the most im- 
portant properties of the ideal gas. The fugacity resp. 
the chemical p oten tial are obtained by a numerical so- 
lution of Eq. (2.6) for given temperature and particle 
number which is exact inasmuch as it does not invoke 
the Thomas-Fermi approximation. The results are then 
analyzed in certain limits below. Fig. l] shows the depen- 
dence of the chemical potential /i on the atom number iV 



for small temperatures. Whereas the solid lines corre- 
spond to the Thomas-Fermi approximation, the three 
other curves wer e obtained numerically by truncating the 
sum in Eq. (2.6) at sufficiently high v. They display a 
step-like variation that becomes increasingly smoother 
for higher temperatures. The step function will appear 
to be broken into smaller steps in the anisotropic oscilla- 
tor case studied below. Here, the steps occur whenever 
a shell is saturated and vp acquires the next higher in- 
teger value, fj, converges to a certain ("plateau") value 
Tiuj(vp + 3/2) in the limit T — > for all TV which do not 
coincide with a "magic" number S a . However, if a shell 
is closed (N = S UF ), fi takes the value \x = hui(vp + 2), 
which is very close to the value of the Thomas-Fermi re- 
sult at S UF (solid line in Fig. 0). As can be shown by 
asymptotic expansion, the two curves intersect approxi- 
mately at N — S a respectively N — (S a + S a -i)/2, i.e. 
at total or half filling of shells. 

This information is displayed in detail in Fig. ^| giving 
the dependence of /i on the temperature around the value 
N = S 7 = 120. At T = Q, the N = 119 curve still 
approaches the previous plateau value (i/hcv = 7 + 3/2, 
whereas the N = 121 curve has to approach the value 
ti/hu} = 8 + 3/2. 

The temperature dependence of /x can be calculated 
analytically in the limits of high and low temperature. 
The high temperature region of Fig. |^ is well described 
by the Sommcrfeld-like formula [Q] 



«"-M'-t(£ 



(2.12) 



where the factor 7r 2 /3 replaces the factor of 7r 2 /12 of the 
usual case of fermions in a box. We note that for high 
temperatures the exact result for N = 119 approaches 
the Sommerfeld approximation for N = 120 (see Fig. |^). 

In the low temperature regime, the variation of /x can 
be analyzed in analogy to the chemical potential of elec- 
trons in an intrinsic semi-conductor. We first consider 
the "magic" numbers N e {S a }. Let N > (T) be the num- 
ber of atoms excited to states above Ep and N < (T) the 
number of unoccupied states ( "holes" ) at or below Ep 



N >( T )= E 7=1 



9v 



^— ' , z 1 eTm(Bfiujy) + 1 

VF f 1 
N<.{T) = ][> [I - z -i exp(/3 ^ ) + 1 

For low temperatures, i.e. for 



(2.13) 
(2.14) 



k B T < E VF+ i - (J, and k B T < n - E F , (2.15) 

the number of "particles" and "holes" can be approxi- 
mated as 



A^<(T) « £ <e -' 3( ' 1 - jEF) , 



(2.16) 
(2.17) 
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where 



C. Specific heat 



£> = ^2 g u e 



-0{E u -E„ F+1 ) 



v=0 



-0{E F -E„) 



(2.18) 
(2.19) 



are essentially Boltzmann sums. On combination of the 
above equations one arrives at 



N > (T)-N < (T)=V > V < e 



-phu 



(2.20) 



From this condition the chemical potential can be deter- 
mined. As for TV s {S a } the Fermi shell is totally filled 
at T = 0, N>(T) must equal N< (T) in that case, and 
one obtains from Eq. fl2.lt -2.2C) the low temperature 
behavior 



M (T) = M^ + 2)-^ln(|^ 



(2.21) 



Thus /j,(0) lies in the middle of the "gap" between Ep 
and Ep + froj, like in an intrinsic semi-conductor where 
the valence band is filled at zero temperature and the 
chemical potential lies in the middle of the bandgap. It 
shows a slow linear decrease with increasing temperature, 
governed by the small factor In (£>/£<). 

If on the other hand N ^ {S a }, one can calculate fi(T) 
from the following approximation. For very low tempera- 
tures the Fermi function is well approximated by n„ = 1 
for v < vp resp. n v — for v > vp , The number of 
occupied states in the Fermi shell AN = N — S UF -i then 
reads approximately 



AN 



9v F 



z 1 exjp(f3hojvp) + 1 



(2.22) 



Assuming that AiV is a constant for very low tempera- 
tures, one can solve Eq. (2.22) for the chemical potential 



lx (T)=?^{vF + -)-kT\n(^-l 



(2.23) 



This expression varies linearly with T for non-vanishing 
AN, with its slope changing sign at AN = g VF /2. If the 
highest shell is less than half filled (AN < g VF /2), fi(T) 
decreases linearly from /i(0) = E VF - if it is more than 
half filled, it increases linearly from the //(0) = E VF . The 
exact result then approaches the Sommerfeld curve. The 
range of validity AT of the linear approximation can be 
roughly determined by equating 



kAT\n(g F - 1) 



(2.24) 



since the maximum deviation from the Sommerfeld ap- 
proximation equals huo/2 at AN = 1 (see Fig.||). For 
vp = 7 one obtains this range as /c^AT '/Tuo < 0.14. For 
larger values of A7V the slope is smaller and the validity 
range may be larger. 



The discontinuity of the chemical potential manifests 
itself most drastically in the specific heat of the gas. It 
is calculated from the total energy 



guTuvv 



z 1 exp(f3hujv) 



C{T) 



dU(T) 
dT ' 



(2.25) 



(2.26) 



The usual Sommerfeld approximation for low tempera- 
tures yields 



C(T) 
Nk 



kT 



hu(<oN) 1 l3 



(2.27) 



whereas the classical high temperature limit equals 
C c ijNk = 3. Here, we determine C(T) for finite N from 
the state sum and compare it to the Sommerfeld approx- 
imation. The results are shown in Fig. ||, where the N 1 ^ 3 
scaling is already included on the ordinate. In the limit 
of ultra-low temperatures the finite-size effects result in 
a deviation from the linear Sommerfeld prediction. For 
higher temperature the s pecific heat C(T) approaches the 
Sommerfeld result (2.27). 

At very low temperatures C (T) remains zero instead of 
increasing linearly. This is consistent with the assump- 
tions leading to Eq. ( 2.23 ) which are confirmed by the cal- 
culation of the state sums. In the ultra-low temperature 
regime no states above the Fermi energy are populated 
due to the energy gap, the total energy does not increase 
and C{T) equals zero. This explanation seems to be cor- 
rect also for the case of closed shells (N = 84, 120, 9880) 
where Eq. ( [2.23 ) does not hold. Note that the range 
where C(T) remains zero is the same as the range of 
validity of the linear approximation for /i(T). 

At intermediate temperatures a strong non-monotonic 
TV-dependence of the specific heat at constant T occurs, 
roughly at those temperatures where the linear approxi- 
mation ceases to be applicable. The origin of this behav- 
ior is revealed in Fig. [|. Each time a shell closure occurs, 
C(N) runs through a maximum. At these points, the 
system can access a new, totally empty shell, at the ex- 
pense of adding the gap energy to the new particles. On 
the contrary, the minima occur half way between succes- 
sive shell closures. In Fig. 0, the total heat capacity is 
plotted versus N without a rescaling. 

Fig. H shows another interesting detail: the two lim- 
iting curves for totally filled shells (N = 84, 120, 9880) 
resp. half filled shells (N = 102, 142, 10270) do not 
depend on up. Up to ksT/hLo m 0.5 the function 
(6A^) 1 / 3 C/(iVA;) does not seem to depend on N explicitly 
for the values considered here, but only on the relative 
filling of the Fermi shell. This is related to the fact that 
the dependence /x(T) shown in Fig. ||repeats itself around 
each value for Pp. 
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D. Density distributions 

Density distributions in traps can be measured quite 
easily. In an isotropic oscillator, one expects radially 
symmetric distributions. The radial wavefunctions u nTt i 
(cf. e.g. ||) are numbered by a radial quantum number 
n r and angular momentum /. The corresponding energy 
is E nrt i = 7Lo>(2(n r -l)+Z + 3/2) so v = 2{n r -l)+l. To 
compute the total density one has to sum up the squared 
wavefunctions weighted correctly with n v 



["/2] + l 



21 + 1 
4tt 1 



Art 



(2.28) 



where l = v - 2(n r - 1). The factor [21 + l)/(47r) is due 
to the summation over all states with m = — I, . . . , I. In 
Fig. [(], p(r) is displayed for different particle numbers and 
temperatures and scaled with the size of the trap ground- 
state a = y/%/(Mu). The zero temperature result from 
the Thomas-Fermi approximation (cf. j^j], see Eq. (3.6), 
A = 1) is also shown. 

For N = 120 (closed shell, vp — 7), one observes a 
central minimum that disappears at N = 142 (half filled 
shell, up = 8). This shell gets totally filled at N = 165 
where p(r) has a maximum at r = 0. The TV-dependence 
of the density at r = is due to the fact that shells with 
odd v do not contribute to p(0) because they are made 
up of odd angular momentum states which all have zero 
density at the origin. The curves for UbT = 0.1?ia; are 
almost indistinguishable from the T = curves. The 
minima and maxima are still visible at kgT = §.25Tilj 
but disappear for temperatures above kgT — hcu. For 
not too high temperatures the density approaches the 
Thomas-Fermi result. Interestingly, p(r) is almost equal 
to this approximation for half filled Fermi shells. 



III. THE ANISOTROPIC CASE 

A. Chemical potential and heat capacity 

Experimentally realized magnetic traps are usually at 
least slightly anisotropic. In this section we therefore 
study a deformed oscillator with a potential 



(3.1) 



i.e. we allow for prolate and oblate ellipsoid iso-energy 
surfaces. Accordingly, the energy eigenvalues are 

E Vr ,u z = hcj(v r + 1 + X(y z + h), (3.2) 

where v r and v z count radial and longitudinal excitations, 
respectively. For given v r there are still g Vr = v T + 1 
degenerate states with different numbers of excitations 
in the two degenerate transversal directions (number of 



partitions of v r = v x + v y ). We use the notation (v r , v z ) 
for these states. 

The TV-dependence of quantities at zero temperature 
turns out to have more features than in the isotropic 
case. For example only in the oblate case (A > 1) the 
notion of shell closures still exists because only then it is 
energetically favorable to fill up a transversal shell with 
degeneracy g Ur before populating a higher longitudinal 
state. However, these new structures occur on a smaller 
scale of particle numbers (due to the smaller degeneracy 
factors) and might be less accessible in experiments with 
a finite uncertainty of the atom number. 

Formulas analogous to Eq. (2.4) can only be given as 
sums. If A < 1, one can count all states up to a certain 
excitation [a r ,a z ) by 



a r -\-[a z A] 



X 



(3.3) 



where [x] denotes the largest integer less than or equal 
to x. Thus, one needs TV = S ari0tz particles to populate 
all states up to E ar , az . For A > 1, there is an analogous 
expression. In general, the exact Fermi energy can only 
be found by searching the lowest state [y ri v z ) with N < 
S„ r>Vr . However, the Thomas- Fermi approximation for 
the Fermi energy is only slightly modified M 



E F = huj(6NX) 1/3 . 



(3.4) 



The Sommerfeld formula Eq. (2.12) for the chemical po- 
tential holds equally in the anisotropic case. 

In the following, we give some of the results for 
anisotropic traps. We concentrate on the case of a heavily 
deformed cigar-shaped trap, A -C 1 (in fact A = 0.076 like 
in [|l0|). We first consider the dependence of the chemi- 
cal potential on the temperature in the low-temperature 
regime. The graph of p(T) for 1000 particles displayed 
in Fig. ^ shows a very intriguing fe ature : it starts lin- 
early at T = as predicted by Eq. (|2.23 ) but then goes 
through a local maximum. The highest occupied s tate 
at T = is (5,22). So g Vr = 6 and with Eq. Q 
AN = 5 5j2 2 - N = 5. The next higher state is (6,9). 
The correction to the linear approximation including this 
state basically shows that it is responsible for the local 
maximum. Indeed, the energy difference between the two 
states corresponds to ksT = 0.157fkoX which is roughly 
at the local minimum (cf. arrow in Fig. [?]). Thus, at the 
maximum of p{T) the next higher level above the Fermi 
level becomes thermally accessible and p decreases. 

As in the isotropic case, the specific heat shows devi- 
ations from the Sommerfeld approximation at tempera- 
tures where the linear approximation for p(T) begins to 
fail. If T is fixed to a value in that region, the graph 
of (6NX)^ 3 C/(Nk) as a function of particle number 
(Fig. ^) exhibits structure on two scales of the particle 
number, considerably more complex than the isotropic 
analogue in Fig. ^. The big jumps take place whenever 
there are enough particles to access a new shell of the 
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transversal oscillator. The arrows denote the values S art o 
for a r = 5, 6, 7. Between two such particle numbers (say 
Sa r ,o, Sa r +i,o) there are 13 major peaks, corresponding 
to 1/A « 13 longitudinal states being filled before the 
next transversal shell can be reached. The finer sub- 
structure (see inset of Fig. ||) can also be explained eas- 
ily: if one starts with the state (a r ,0) the next state is 
(0, [ov/A] + 1) followed by (1, [(a r - 1)/A] + 1), etc. Con- 
sequently, there should be a r maxima before (a r , 1) is 
reached. In case 1/ A is integer this substructure disap- 
pears. 

B. Density distributions 

In order to calculate the density distribution we make 
use of the transverse symmetry and obtain 

oo K/2] 
p(r,z) = 2 ^ n v r ,vS T ) / , \Un r ,Vr-2n r {r) ■ 1pn z (z)\ 2 , 
Vr,i>'z=0 n r —0 

(3.5) 

where u nr ,u r -2n T (i r ) is the radial wavefunction of a 2D 
harmonic oscillator (r 2 = x 2 + y 2 ) with magnetic quan- 
tum number \m\ = v r — 2n r and ipn z { z ) is the wave- 
function of the one-dimensional harmonic oscillator. The 
overall factor 2 allows for the twofold degeneracy of a 
state with given v z ,N r , \tn\. n Vr ,u z {T) is the Fermi-Dirac 
occupation number analogous to Eq. ( |2.5|) . The numer- 
ical results should again be compared to the Thomas- 
Fermi approximation at T = 

NX 8 / rW\ 3/2 
P{rt * ) = ^[ 1 -^—) ' (3 - 6) 

where Rf = (48A^A) 1 / 6 cr is the so called "Fermi radius" 
of the density distribution. 

We restrict ourself to a plot of p(r, 0) resp. p(0,z) 
(cf. Fig. ||). In transversal direction the density at zero 
temperature shows only very slight deviations from the 
Thomas- Fermi result. In contrast, the density in lon- 
gitudinal direction shown on the right exhibits oscilla- 
tions around this approximation. These oscillations are 
mainly due to a step- like behavior of the occupation num- 
ber N Vx of the oscillator states in longitudinal direction, 
which is due to the filling of new transversal states. The 
steps have a width of 1/A. In addition, N Vm decreases 
much slower (average derivative oc A) than the occu- 
pation number N Vx for the transversal states (average 
derivative oc 1/A, here we consider p(x, y = z — 0)). As 
a result, the longitudinal density profile p(0, z) receives a 
bigger contribution from higher states than the transver- 
sal one, so one can still see the various maxima of high 
oscillator states in the first case but not in the latter. 

We finally note that for oblate traps the effects are 
basically the same as in prolate traps. The behavior of 



the density profiles is interchanged, the abovementioned 
oscillations occur in the transversal density profile and 
disappear in the longitudinal direction. 

IV. CONCLUSION 

Our calculations have shown that quantum statistical 
effects on the easily accessible observables of a trapped 
Fermi gas are restricted to the regime of rather small 
atom number, e.g. below N — 1000. For larger atom 
numbers quantum statistical effects can be more easily 
understood in terms of a local density approximation 
0,0. Our work has been carried out using the grand 
canonical ensemble. It is known from the ideal Bose gas 
that the grand canonical and the canonical ensemble give 
differing predictions ]§| . However the deviations are often 
smaller than the difference between the interacting and 
the interaction-free case. Problems like artificial fluctu- 
ations in a grand canonical ensemble of bosons will not 
occur for fcrmions. For that reason we don't dwell on a 
comparison of the different ensembles here. 

The present study of an ideal Fermi gas showed some 
remarkable effects of the shell structure in the harmonic 
potential visible e.g. in Fig. [|. For a real Fermi gas the 
atom-atom interactions will introduce an additional de- 
pendence of the chemical potential on the particle num- 
ber which will smoothen out the steps in Fig. [j]. A quan- 
titative prediction of this effect depends on the relative 
magnitude of atom-atom interactions and the experimen- 
tally controllable energy gap huj. As mentioned in the in- 
troduction, an effective suppression of the interactions in 
the ultracold regime can be attributed to the suppression 
of s-wave scattering for Fermions in identical substates. 
Therefore the ideal gas behavior might still be visible 
provided the harmonic potential is steep enough. 

BCS-like behavior as discussed in [|||| requires two 
spin states to be trapped. The effects of the harmonic 
potential in that situation have been allowed for in local 
density approximation in ||, i.e. for large atom numbers. 
Because of the sensitivity of the BCS transition to the 
difference of the atom numbers in both spin states, it 
however remains a challenge to observe a BCS transition 
experimentally. By contrast shell effects as discussed in 
this paper should be visible in a suitable range of small 
atom numbers and sufficiently large trap frequencies. 
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FIG. 1. /x(TV) shows a step-l ike b ehavior following the con- 
tinuous approximation in Eq. ( |2,ll| ). The dashed and dotted 
curves are displaced vertically by —1 resp. —2. 
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FIG. 2. This figure shows /x(T) for the isotropic trap. The 
Sommerfeld approximation is for TV = SV = 120 but agrees 
very well with the numerical curve for TV = 119. This occurs 
also for other values of TV. 
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FIG. 3. C(T) appro priate ly scaled. The linear curve is the 
Sommerfeld result Eq. (2.27). The deviation at large fcjl '/hcj 
from the linear behavior occurs only for small TV, because for 
high temperatures (6TV) 1/3 C(T)/(Nk) -> 3(6TV) 1/3 due to 
the equipartition theorem. 
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FIG. 4. Specific heat as a function of the particle num- 
ber at different temperatures and scaled with (67V) 1 / 3 . 
r = kBT/hui denotes the temperature in units of the level 
spacing. The arrows point to TV = S v for v = 1, . . . , 11 
(N = 4, 10, 20, 35, 56, 84, 120, 165, 220, 286, 364). 
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FIG. 5. Total specific heat as a function of the particle 
number at different temperatures. Again, r = fcsT /hco. 
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FIG. 6. Spatial density p in a isotropic trap as a function 
of the distance to the trap center for different numbers of 
fermions. The scaling parameters a — ^jTi/(Mu>) is the width 
of the ground state. The unbroken lines denote ksT = Q.lhu), 
the dotted fesT = 0.25?ki; and the dashed ones are obtained 
from the Thomas-Fermi approximation at T = 0. 
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FIG. 7. fJ,(T) for a cigar-shaped trap with A = 0.076 and 
TV = 1000. The arrow denotes the energy difference between 
the Fermi level and the next higher state. The inset displays 
the linear and next higher approximation taking into account 
the states directly below and above the Fermi level. 
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FIG. 8. Specific heat for a cigar-shaped trap as a func- 
tion of TV at k B T/%Lj = 0.044. The arrows point to 
Sa r ,o (ov = 5,6,7). The inset shows the number range 
750 < TV < 850 in more detail. 
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FIG. 9. Density p for a cigar-shaped trap. The left part 
shows the density in transversal direction for 2 = 0, the graph 
on the right is the density in longitudinal direction on the 
symmetry axis (x = y — 0). The unbroken lines are nu- 
merical results for T = 0, the dashed lines come from the 
Thomas-Fermi approximation. Note that the "Fermi radius" 
in longitudinal direction is about 1/A ~ 13 times larger than 
the transversal one. 
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